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Abstract 

We provide a new method for constructing equiangular tight frames (ETFs). The con- 
struction is valid in both the real and complex settings, and shows that many of the few 
previously-known examples of ETFs are but the first representatives of infinite families of 
such frames. It provides great freedom in terms of the frame's size and redundancy. This 
method also explicitly constructs the frame vectors in their native domain, as opposed to 
implicitly defining them via their Gram matrix. Moreover, in this domain, the frame vec- 
tors are very sparse. The construction is extremely simple: a tensor-like combination of a 
Steiner system and a regular simplex. This simplicity permits us to resolve an open question 
regarding ETFs and the restricted isometry property (RIP): we show that the RIP behavior 
of some ETFs is unfortunately no better than their coherence indicates. 
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1. Introduction 

Let F = {f n }n =l be a finite sequence of vectors in a real or complex M-dimensional 
Hilbert space M M . The corresponding frame operator is FF* = 52 n=1 fnfni where /* denotes 
the linear functional that maps a given / G Mm to the scalar (/, f n ). The sequence F is said 
to be a tight frame if there exists A > such that FF* = Al. Meanwhile, F is equiangular 
if || / n || = 1 for all n and if there exists a > such that |(/ n , f n >) \ = a for all n ^ n' . This 
paper concerns equiangular tight frames (ETFs); writing FasanMxiV matrix, we need the 
rows of F to be orthogonal and have constant norm, the columns of F to be unit norm, and 
the inner products of distinct columns of F to have constant modulus. As detailed below, 
such frames are useful in applications, but up to this point, they have proven notoriously 
difficult to construct. 

In this article, we provide a new method for constructing ETFs. The construction is 
valid in both the real and complex settings, and shows that many of the few previously- 
known examples of ETFs are but the first representatives of infinite families of such frames. 
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This construction technique also permits great freedom in selecting M and N, just shy of 
letting one choose the exact size and redundancy of their liking. This method also explicitly 
constructs the frame vectors in their native domain as opposed to the usual method of 
implicitly defining them with their Gram matrix F*F. Moreover, in this domain, the frame 
vectors can be chosen to be very sparse. The construction is extremely simple: a tensor-like 
combination of a Steiner system and a regular simplex. This simplicity permits us to resolve 
an open question regarding ETFs and the restricted isometry property (RIP): we show that 
the RIP behavior of some ETFs is unfortunately no better than their worst- case- coherence 
bounds indicate. 

Equiangular lines have long been a subject of interest 19|. Recent work on the matter of 
ETFs was spurred on by communications-theory- inspired results 0, 16, 27] that show that 
the linear encoders provided by such frames are optimally robust against channel erasures. 
In the real setting, the existence of an ETF of a given size is equivalent to the existence 



of a strongly regular graph with certain corresponding parameters [la, |24|. Such graphs 
have a rich history and remain an active topic of research [f|; the specific ETFs that arise 



from particular graphs are detailed in 30]. Some of this theory gen eralizes to the complex- 



variaUe setting in the guise ot compfa Seidel matrices g &U Many approaches to 
constructing ETFs have focused on the special case in which every entry of F is a root of 



unity [17|, |20j, |26|, |28|, |3jj . Other approaches are given in [n|, |25|, |29| . In the complex setting, 
much attention has focused on the maximal case of M 2 vectors in M M 0, 18, 21, 23 . 

A version of the ETF construction method we present here was previously employed by 
Seidel in Theorem 12.1 of 24] to prove the existence of certain strongly regular graphs. In 
the context of that result, our contributions are: (i) the realization that when Seidel's block 
design arises from a particular type of Steiner system, the resulting strongly regular graph 
indeed corresponds to a real ETF; (ii) noting that in this case, the graph theory may be 
completely bypassed, as the idea itself directly produces the requisite frame F; and (iii) 
having bypassed the graph theory, realizing that this construction immediately generalizes 
to the complex-variable setting if Seidel's requisite Hadamard matrix is permitted to become 



12 



complex. These realizations permit us to exploit the vast literature on Steiner systems 
to construct several new infinite families of ETFs, in both the real and complex settings 
Moreover, these ETFs are extremely sparse in their native space; sparse tight frames have 
recently become a subject of interest in their own right 
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In fact, these ETFs are simple enough so as to permit a rigorous investigation of their 
potential as RIP matrices, which are currently in demand due to their applicability in 
compressed sensing 0, @] . As discussed below, ETFs are the optimal matrices with respect 
to a very coarse estimate — worst-case coherence — on a matrix's RIP bounds. Our hope was 
that all ETFs, having such high degrees of symmetry, might possess other hidden properties 
that, when properly exploited, yield even better bounds than those given by coherence-based 
estimates. Unfortunately, our newly-discovered ETF constructions dash these hopes: for at 
least some ETFs, worst-case coherence, a la Gershgorin circles, does indeed provide a very 
good estimate on RIP bounds. With respect to RIP, these ETFs perform no better than a 
myriad of previously discovered deterministic constructions of RIP matrices, such as those 



given m 



131. 



In the next section, we provide our main result, namely Theorem HJ which shows how 
certain Steiner systems may be combined with regular simplices to produce ETFs. In the 
third section, we discuss each of the known infinite families of such Steiner systems, and 
compute the corresponding infinite families of ETFs they generate. We further provide 
some necessary and asymptotically sufficient conditions, namely Theorem [2J to aid in the 
quest for discovering other examples of such frames that lie outside of the known infinite 
families. In Section 4, we discuss the possible RIP behavior of ETFs in general, and show 
that the performance of our Steiner ETFs is indeed no better than that guaranteed by 
coherence-based estimates. 



2. Steiner equiangular tight frames 

In this section, we provide new constructions of infinite families of ETFs, namely M x N 
matrices F = [f\... Jn] which have orthogonal rows of constant squared-norm A and unit 
norm columns whose inner products have constant modulus a: we want FF* = Al while 
the diagonal entries of F*F are 1 and the off-diagonal entries are a in modulus. For a fixed 
M and N, there is no ambiguity [27| as to the values of A and a. Indeed, noting that since 
N = J2n=i ll/nll 2 = Tr(F*F) = Tr(FF*) = MA, we have A = §; moreover, since 

N 



N + N(N-l)a 2 = \(f n Jn')\ 2 = TT[(F*F) 2 } = Tr[(FF*) 2 ] = MA 



2 _ N± 
M ' 



then a 2 = jgfejrn • Conversely, if one can design an N x N self-adjoint, positive semidefinite 
Gram matrix G of rank M whose diagonal entries are one and whose off-diagonal entries are 
an m[ 7j M -\ i i n squared- modulus, then one can then factor G as F*F, where F is an M x N 



i) 



ETF [27|. This fact has led many to attempt to construct ETFs, not by constructing F 
directly, but rather, by constructing G. This Gram representation of an ETF has the addi- 
tional benefit of being invariant with respect to rotations of the frame elements themselves, 
and, in the real- variable case, is closely-related to the incidence matrix of the corresponding 



strongly regular graph [16] . Moreover, whenever G = F*F and FF* = Al, then the columns 
of G are, in fact, a scalar multiple of an isometric embedding of the columns of F. That is, 
the columns of the Gram matrix of a tight frame are but a high-dimensional representation 
of the frame elements themselves. 

There is a drawback, however, to working with Gram representations: one does not pro- 
duce the frame vectors in their native M-dimensional space, the domain in which they are 
usually needed for communications applications. And though factoring G is straightforward — 
one may, for instance, apply the Gram-Schmidt algorithm to the columns of G — it will 
produce the / n 's with respect to some arbitrarily chosen basis for one that may not 
be optimal. Indeed, such a process ignores the question of whether or not there is a basis 
for Mm that makes the frame elements sparse. For this reason, in this paper, we avoid the 
Gram representation and construct F directly. The key idea is to design the ETFs in blocks, 
specifically those arising from a particular type of combinatorial block design. 
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Steiner systems and block designs have been studied for over a century; the background 
facts presented here on these topics are taken from [l|, In short, a {y,b,r,k,X)-block 
design is a v -element set V along with a collection B of b fc-element subsets of V, dubbed 
blocks, that have the property that any element of V lies in exactly r blocks and that any 
2-element subset of V is contained in exactly A blocks. The corresponding incidence matrix 
is a v x b matrix A that has a one in a given entry if that block contains the corresponding 
point, and is otherwise zero; in this paper, it is more convenient for us to work with the 
b x v transpose A T of this incidence matrix. Our particular construction of ETFs involves 
a special class of block designs known as (2, k, v)-Steiner systems. These have the property 
that any 2-element subset of V is contained in exactly one block, that is, A = 1. With 
respect to our purposes, the crucial facts are the following: 

The transpose A T of the {0, l}-incidence matrix A of a (2, fc, t> )-Steiner system: 

(i) is of size x v, 

(ii) has k ones in each row, 

(iii) has jrjr ones in each column, and 

(iv) has the property that any two of its columns have a dot product of one. 

The first three facts follow immediately from solving for b = and r = |5y, usrn g the 

well-known relations vr = bk and r(k — 1) = X(v — 1). Meanwhile, (iv) comes from the fact 
that A = 1: each column of A T corresponds to an element of the set, and the dot product 
of any two columns computes the number of blocks that contains the corresponding pair of 
points. This in hand, we present our main result; here, the density of a matrix is the ratio 
of the number of nonzero entries of that matrix to the entire number of its entries: 

Theorem 1. Every (2, k,v)- Steiner system generates an equiangular tight frame consisting 
of N = + vectors in M = - dimensional space with redundancy |j = k(l + |5j0 
and density J = ( M( ^~ 1 M) )i 

Moreover, if there exists a real Hadamard matrix of size 1 + then such frames are real. 
Specifically, a x v(l + ETF matrix F may be constructed as follows: 

1. Let A T be the x v transpose of the adjacency matrix of a (2, k, v)-Steiner system. 

2. For each j = 1, . . . , v , let Hj be any (1 + f^f) x (1 + |5y) matrix that has orthogonal 
rows and unimodular entries, such as a possibly complex Hadamard matrix. 

3. For each j = 1, . . . ,v , let Fj be the jjfe^ x (1 + matrix obtained from the jth 
column of A T by replacing each of the one-valued entries with a distinct row of Hj, 
and every zero-valued entry with a row of zeros. 

4. Concatenate and rescale the Fj 's to form F = (^zj)^[Fi ■ ■ ■ F v ]. 

We refer to the ETFs produced by Theorem [T] as (2, k, v) -Steiner ETFs. In essence, the 
idea of the construction is to realize that the nonzero rows of any particular Fj form a regu- 
lar simplex in |5i"dimensional space; these vectors are automatically equiangular amongst 

4 



themselves; by requiring the entries of these simplices to be unimodular, and requiring that 
distinct blocks have only one entry of mutual support, one can further control the inner 
products of vectors arising from distinct blocks. This idea is best understood by consid- 
ering a simple example, such as the ETF that arises from a (2, 2, 4)-Steiner system whose 
transposed incidence matrix is: 



A 1 
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+ 
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+ 
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+ 



One can immediately verify that A T corresponds to a block design: there is a set V of v — 4 
elements, each corresponding to a column of A T ; there is also a collection B of h = 6 subsets 
of V, each corresponding to a row of A T ; every row contains k = 2 elements; every column 
contains r = 3 elements; any given pair of elements is contained in exactly one row, that is, 
A = 1, a fact which is equivalent to having the dot product of any two distinct columns of 
A T being one. To form an ETF, for each of the four columns of A T we must choose a 4 x 4 
matrix H with unimodular entries and orthogonal rows; the size of H is always one more 
than the number r of ones in a given column of A T . Though in principle one may choose 
a different H for each column, we choose them all to be the same, namely the Hadamard 
matrix: 

+ + +" 
H " + " 



To form the ETF, for each column of A T we replace each of its 1-valued entries with a 
distinct row of H . Again, though in principle one may choose a different sequence of rows 
of H for each column, we simply decide to use the second, third and fourth rows, in that 
order. The result is a real ETF of N = 16 elements of dimension M = 6: 



F= —= 



1 



+ - 
+ + 



+ 



+ 



+ + 
+ - 



+ 



-- + + -- 



+ - + 



+ + - - 



One can immediately verify that the rows of F are orthogonal and have constant norm, 
implying F is indeed a tight frame. One can also easily see that the inner products of two 
columns from the same block are — |, while the inner products of columns from distinct 
blocks are ±|. Theorem [T] states that this behavior holds in general for any appropriate 
choice of A T and H; its formal proof is as follows. 
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Proof. To verify F is a tight frame, note that the inner product of any two distinct rows of 
F is zero, as they are the sum of the inner products of the corresponding rows of the Fj's 
over all j = 1, . . . ,v; for any j, these shorter inner products are necessarily zero, as they 
either correspond to inner products of distinct rows of Hj or to inner products with zero 
vectors. Moreover, the rows of F have constant norm: as noted in (ii) above, each row of 
A T contains k ones; since each Hj has unimodular entries, the squared-norm of any row of 
F is the squared-scaling factor times a sum of k(l + rz|) ones, which, as is necessary 
for any unit norm tight frame, equals the redundancy |j = k(l + frr)- 

Having that F is tight, we show F is also equiangular. We first note that the columns 
of F have unit norm: the squared-norm of any column of F is times the squared- 
norm of a column of one of the Fj's; since the entries of Hj are unimodular and (iii) above 
gives that each column of A T contains |5y ones, the squared-norm of any column of F is 
(^zf)(fzi)l = 1' as claimed. Moreover, the inner products of any two distinct columns of F 
has constant modulus. Indeed, the fact (iv) that any two distinct columns of A T have but 
a single entry of mutual support implies the same is true for columns of F that arise from 
distinct Fj blocks, implying the inner product of such columns is times the product 
of two unimodular numbers. That is, the squared-magnitude of the inner products of two 
columns that arise from distinct blocks is j^r^n = (^E~[) 2 , as needed. Meanwhile, the same 
holds true for columns that arise from the same block Fj. To see this, note that since Hj is 
a scalar multiple of a unitary matrix, its columns are orthogonal. Moreover, Fj contains all 
but one of the Hj's rows, namely one for each of the 1-valued entries of A T , a la (iii). Thus, 
the inner products of the portions of Hj that lie in Fj are their entire inner product of zero, 
less the contribution from the left-over entries. Overall, the inner product of two columns 
of F that arise from the same Fj block is times the negated product of one entry of Hj 
and the conjugate of another; since Hj is unimodular, we have that the squared-magnitude 
of such inner products is jjffej^n = (f~[) 2 ; as needed. 

Thus F is an ETF. Moreover, as noted above, its redundancy is |r = k(l + ^5? )■ All that 
remains to verify is its density: as the entries of each Hj are all nonzero, the proportion of 
F's nonzero entries is the same as that of the incidence matrix A, which is clearly having k 
ones in each v -dimensional row. Moreover, substituting N — v(l + |^y) and M = jjfe^ into 
the quantity m ^~ 1 m) reveals it to be ^, and so the density can be alternatively expressed 
as ( M(N-M) ) * > as claimed. □ 

In the next section, we apply Theorem [TJ to produce several infinite families of Steiner 
ETFs. Before doing so, however, we pause to remark on the redundancy and sparsity of such 
frames. In particular, note that since the parameters k and v of the requisite Steiner system 
always satisfy 2 < k < v, then the redundancy k(l + ^5r) °f Steiner ETFs is always between 
k and 2k; the redundancy is therefore on the order of k, and is always strictly greater than 
2. If a low-redundancy ETF is desired, one can always take the Naimark complement [5| 
of an ETF of N elements in M-dimensional space to produce a new ETF of N elements in 
(N — M)-dimensional space; though the complement process does not preserve sparsity, it 
nevertheless transforms any Steiner ETF into a new ETF whose redundancy is strictly less 
than 2. However, such a loss of sparsity should not to be taken lightly. Indeed, the low 
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density of Steiner ETFs gives them a large computational advantage over their non-sparse 
brethren. 

To clarify, the most common operation in frame-theoretic applications is the evaluation 
of the analysis operator F* on a given / e M.m- For a non-sparse F, this act of com- 
puting F*f requires O(MN) operations; for a frame F of density D, this cost is reduced 
to O(DMN). Indeed, using the explicit value of D = ( ^(^-m) ) 5 §i ven m Theorem [T] as 
well as the aforementioned fact that the redundancy of such frames necessarily satisfies 

> 2, we see that the cost of evaluating F*f when F is a Steiner ETF is on the order 
°f ( ^jv-M^ )^^" < (2M)^iV operations, a dramatic cost savings when M is large. Further 
efficiency is gained when F is real, as its nonzero elements are but a fixed scaling factor times 
the entries of a real Hadamard matrix, implying F*f can be evaluated using only additions 
and subtractions. The fact that every entry of F is either or ±1 further makes real Steiner 
ETFs potentially useful for applications that require binary measurements, such as design 
of experiments. 



3. Examples of Steiner equiangular tight frames 

In this section, we apply Theorem [T] to produce several infinite families of Steiner ETFs. 
When designing frames for real-world applications, three considerations reign supreme: size, 
redundancy and sparsity. As noted above, every Steiner ETF is very sparse, a serious 
computational advantage in high- dimensional signal processing. Moreover, some of these 
infinite families, such as those arising from finite affine and projective geometries, provide 
one great flexibility in choosing the ETF's size and redundancy. Indeed, these constructions 
provide the first known guarantee that for a given application, one is always able to find 
ETFs whose frame elements lie in a space whose dimension matches, up to an order of 
magnitude, that of one's desired class of signals, while simultaneously permitting one to 
have an almost arbitrary fixed level of redundancy, a handy weapon in the fight against 
noise. To be clear, recall that the redundancy of a Steiner ETF is always strictly greater 



than 2. Moreover, as general bounds on the maximal number of equiangular lines [19 



require that any ETF satisfy iV < MM±D j n rea j spaces and N < M 2 in complex ones, the 
redundancy of an ETF is never truly arbitrary. Nevertheless, if one does prescribe a given 
desired level of redundancy in advance, the Steiner method can produce arbitrarily large 
ETFs whose redundancy is approximately the prime power nearest to the sought-after level. 

3.1. Infinite families of Steiner equiangular tight frames 

We now detail eight infinite families of ETFs, each generated by applying Theorem [T] to 
one of the eight completely understood infinite families of (2, k, v )-Steiner systems. Tabled] 
summarizes the most important features of each family, while Table [5] gives the first few 
examples of each type, summarizing those that lie in 100 dimensions or less. 

3.1.1. All two-element blocks: (2, 2, v) -Steiner ETFs for any v > 2. 

The first infinite family of Steiner systems is so simple that it is usually not discussed 
in the design-theory literature. For any v > 2, let V be a ^-element set, and let B be the 
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v(v— 1) 



blocks, each of which 



collection of all 2-element subsets of V. Clearly, we have b — ,, 
contains k = 2 elements; each point is contained in r = v — 1 blocks, and each pair of points 
is indeed contained in but a single block, that is, A = 1. 

By Theorem [TJ the ETFs arising from these (2, 2, v )-Steiner systems consist of N = 



v 



V 



vectors in M 



v(v—l) v(v— 1) 



dimensional space. Though these frames 



fc-lV " »^^" ^« fc(fc-l) 2 

can become arbitrarily large, they do not provide any freedom with respect to redundancy: 
= 2^- is essentially 2. These frames have density - = -. Moreover, these ETFs can 
be real-valued if there exists a real Hadamard matrix of size 1 + ^Ej — v - I n particular, it 
suffices to have v to be a power of 2; should the Hadamard conjecture prove true, it would 
suffice to have v divisible by 4. 

One example of such an ETF with v — 4 was given in the previous section. For an- 
other, consider v = 3. The b x v transposed incidence matrix A T is 3 x 3, with each row 
corresponding to a given 2-element subset of {0, 1, 2}: 



+ + 
+ + 
+ + 



To form the corresponding 3x9 ETF F, we need a 3 x 3 unimodular matrix with orthogonal 
rows, such as a DFT; letting uj = e 27rl//3 , we can take 



H 



1 1 
1 
1 



CO 



1 

u 2 



To form F, in each column of A T , we replace each 1-valued entry with a distinct row of H. 
Always choosing the second and third rows yields an ETF of 9 elements in C 3 : 



UJ 



UJ UJ 



1 UJ 



UJ 



UJ 
UJ 



UJ 
UJ 2 



This is the only known instance of when the Steiner-based construction of Theorem [TJ pro- 
duces a maximal ETF, namely one that has iV = M 2 . 

3.1.2. Steiner triple systems: (2, 3, v)-Steiner ETFs for any v = 1, 3 mod 6. 

Steiner triple systems, namely (2, 3, t>)-Steiner systems, have been a subject of interest 

Each 
- 2=1 
3 



for over a century, and are known to exist precisely when v = 1,3 mod 6 [12 
of the b = ^""^ blocks contains k = 3 points, while each point is contained in 



vectors in 



blocks. The corresponding ETFs produced by Theorem [TJ consist of ^±±1 , ( 

' As with ETFs stemming from 2-element 



v-1 



is always 



dimensional space. The density of such frames is -. 

blocks, Steiner triple systems offer little freedom in terms of redundancy: , , 
approximately 3. Such ETFs can be real if there exists a real Hadamard matrix of size 

The Fano plane is a famous example of such a design. The simplest example of a finite 
projective geometry, it consists of v = 7 points and 6 = 7 lines, any two of which intersect in 
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exactly one point. Each line consists of k = 3 points, and each point is contained in r 
lines: 



A 1 
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+ 








+ 






+ + 






+ 
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+ 




+ 




+ 


+ 






+ 




+ 




+ 






+ 


+ 




+ 



+ 



Choosing H to be the standard 4x4 Hadamard matrix used in the previous section results 
in a real ETF of 28 elements in 7-dimensional space; F is the scaling factor times: 



3.1.3. Four element blocks: (2, A, v) -Steiner ETFs for any v = 1,4 mod 12. 

It is known that (2, 4, v )-Steiner systems exist precisely when v = 1, 4 mod 12 Q. Con- 
tinuing the trend of the previous two families, these ETFs can vary in size but not in 



redundancy: they consist of 



v(v+2) 
3 



vectors in 



v(v— 1) 



12 



dimensional space, having redundancy 



4^| and density -. Interestingly, such frames can never be real: with the exception of the 
trivial lxl and 2x2 cases, the dimensions of all real Hadamard matrices are divisible by 



4; since v = 1,4 mod 12, the requisite matrices H here are of size 



1,2 



mod 4. 



3.1.4- Five element blocks: (2, 5, v)-Steiner ETFs for any v = 1,5 mod 20. 

It is known that (2, 5, v )-Steiner systems exist precisely when v = 1,5 mod 20 



The 



corresponding ETFs consist of Hlld^i 

and density 
of size ^ 



vectors in ( 2Q -dimensional space, having redundancy 
-. Such frames can be real whenever there exists a real Hadamard matrix 

V 

In particular, letting v = 45, we see that there exists a real Steiner ETF of 540 
vectors in 99-dimensional space, a fact not obtained from any other known infinite family. 

3.1.5. Affine geometries: (2, q, q n ) -Steiner ETFs for any prime power q, n > 2. 

At this point, the constructions depart from those previously considered, allowing both 
k and v to vary. In particular, using techniques from finite geometry, one can show that 
for any prime power q and any n > 2, there exists a (2, k, u)-Steiner system with k = q 



and v = q n 12 



1 + 



_i- 

-1 



vectors in q 



n-Li q n -l - 



The corresponding ETFs consist of q r \- ■ • 

dimensional space. Like the preceding four classes of Steiner ETFs, these frames can grow 
arbitrarily large: fixing any prime power q, one may manipulate n to produce ETFs of 
varying orders of magnitude. However, unlike the four preceding classes, these affine Steiner 
ETFs also provide great flexibility in choosing one's redundancy. That is, they provide the 
ability to pick M and N somewhat independently. Indeed, the redundancy of such frames 
q(l + -^pr) is essentially q, which may be an arbitrary prime power. Moreover, as these 



frames grow large, they also become increasingly sparse: their density is 
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Because of 



their high sparsity and flexibility with regards to size and redundancy, these frames, along 
with their projective geometry-based cousins detailed below, are perhaps the best known 
candidates for use in ETF-based applications. Such ETFs can be real if there exists a real 
Hadamard matrix of size 1 + ^-£j-, such as whenever q = 2, or when q = 5 and n = 3. 

3.1.6. Projective geometries: (2, q + 1, 9 Steiner ETFs for any prime power q, n > 2. 
With finite geometry, one can show that for any prime power q and any n > 2, there 

exists a (2, fc,u)-Steiner system with k = q + 1 and v = 9 _~ 1 [l2f| . Qualitatively speaking, 
the ETFs these projective geometries generate share much in common with their affmely- 
generated cousins, possessing very high sparsity and great flexibility with respect to size and 
redundancy. The technical details are as follows: they consist of 9 Jf 1 ^- + 2 -£j0 vec tors 
i n -dimensional space, with density ^l+rzi an d redundancy (q + 1)(1 + ^rrj)- 

These frames can be real if there exists a real Hadamard matrix of size 1 + 2 - z r; note this 

q—l ' 

restriction is identical to that for ETFs generated by afnne geometries for the same q and n, 
implying that real Steiner ETFs generated by finite geometries always come in pairs, such 
as the 6 x 16 and 7 x 28 ETFs generated when q = 2, n = 2, and the 28 x 64 and 35 x 120 
ETFs generated when q = 2, n = 3. 

3.1.7. Unitals: (2, q + 1, q 3 + 1) -Steiner ETFs for any prime power q. 

For any prime power q, one can show that there exists a (2, k, v )-Steiner system with 
k = q + 1 and v = q 3 + 1 [12]. Though one may pick a redundancy of one's liking, such a 
choice confines one to ETFs of a given size: they consist of (q 2 + 1) (q 3 + 1) vectors in 9 ^qj-^ - 
dimensional space, having redundancy (q + 1)(1 + 4?) and density These ETFs can 

never be real: the requisite Hadamard matrices are of size q 2 + 1 which is never divisible by 
4 since and 1 are the only squares in Z4. 



5.1.5. Denmston designs: (2, 2 r , 2 r+s + 2 r - 2 s ) -Steiner ETFs for any 2 < r < s. 

For any 2 < r < s, one can show that there exists a (2, k, v )-Steiner system with k = 2 r 
and i> = 2 r+s + 2 r — 2 s [l2[. By manipulating r and s, one can independently determine 
the order of magnitude of one's redundancy and size, respectively: the corresponding ETFs 
consist of (2 s + 2)(2 r+s + 2 r — 2 s ) vectors in ( 2 +1 )( 2 + 2 ~ 2 ) -dimensional space, having redun- 
dancy 2 r |Ji| and density 2r+a _^, r _ 2s . As such, this family has some qualitative similarities 
to the familes of ETFs produced by affine and projective geometries. However, unlike those 
families, the ETFs produced by Denniston designs can never be real: the requisite Hadamard 
matrices are of size 2 s + 2, which is never divisible by 4. 

3.2. Necessary and sufficient conditions on the existence of Steiner ETFs. 

(2, k, v )-Steiner systems have been actively studied for over a century, with many cele- 
brated results. Nevertheless, much about these systems is still unknown. In this subsection, 
we discuss some known partial characterizations of the Steiner systems which lie outside of 
the eight families we have already discussed, as well as what these results tell us about the 
existence of certain ETFs. To begin, recall that, for a given k and v, if a (2, k, u)-Steiner 
system exists, then the number r of blocks that contain a given point is necessarily t^Tj 
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Name 



M 



N 



v (v — 1) 
2 

v(v — 1) 



Redundancy 



Real? 



Restrictions 



2- blocks 

3- blocks 

4- blocks 

5- blocks 
Affine 
Projective 
Unitals 
Denniston 



V 

v(v + l) 



6 


2 


v(v — 1) 


l>(u+2) 


12 


3 


— 1) 


v(v+3) 


20 


4 




+ ■ 


( g "-i)( 9 "+ i -i) 





i 



I v+2 



■v+3 



q(l 



v + 1 
2 

Never 

v+3 
4 

1 + 



q-l 



8+1 

(2' + l)(2 T '+ s +2'--2') 

2!: 



q-l v g-1 ' 

(g 2 + l)(q 3 + 1) 
(2 s + 2)(2 r+s + 2 r - 2 s ) 



(g + i)(i + ^) 



2 S +1 



9- 

Never 
Never 



None 

d = 1, 3 mod 6 

v = 1,4 mod 12 

v = 1, 5 mod 20 

g a prime power, n > 2 

g a prime power, n > 2 

ij a prime power 

2 < r < s 



Table 1: Eight infinite families of Steiner ETFs, each arising from a corresponding known infinite family of 
(2, k, w)-Steiner designs. Each family permits both M and N to grow very large, but only a few families — 
affine, projective and Denniston — give one the freedom to simultaneously control the proportion between M 
and N, namely the redundancy ^ of the ETF. The column denoted "Real?" indicates the size for which a 
real Hadamard matrix must exist in order for the resulting ETF to be real; it suffices to have this size be a 
power of 2; if the Hadamard conjecture is true, it would suffice for this number to be divisible by 4. 



while the total number of blocks b is ^fezjl- ^ s such, in order for a (2, k, v )-Steiner system 
to exist, it is necessary for (k, v) to be admissible, that is, to have the property that |5y and 

V(V— 1) ■ ; 

yj—fj are integers. 

However, this property is not sufficient for existence: it is known that a (2, 6, 16)-Steiner 
system does not exist [l[ despite the fact that |5i — 3 an d kQc-i) = ^ I n ^ ac ^' letting v be 
either 16, 21, 36, or 46 results in an admissible pair with k = 6, despite the fact that none of 
the corresponding Steiner systems exist; there are twenty- nine additional values of v which 
form an admissible pair with k = 6 and for which the existence of a corresponding Steiner 
system remains an open problem [l| . Similar nastiness arises with k > 7. The good news is 
that admissibility, though not sufficient for existence, is, in fact, asymptotically sufficient: 
for any fixed k, there exists a corresponding admissible index Vo(k) for which for all v > Vo(k) 
such that and jfe^ are integers, a (2, k, u)-Steiner system indeed exists jlj. Moreover, 
explicit values of vo(k) are known for small k: t>o(6) = 801, v (7) = 2605, t>o(8) = 3753, 
t>o(9) = 16497. We now detail the ramifications of these design-theoretic results on frame 
theory: 

Theorem 2. // an N -element Steiner equiangular tight frame exists for an M -dimensional 
space, then letting a = ( m(n-i) ) ^ > ^ e corres P° n ding block design has parameters: 

»• - Na u — M, r = K k- N 



l+a' ' a' M(l+a)^ 

In particular, if such a frame exists, then these expressions for v, k and r are necessarily 
integers. 

Conversely, for any fixed k > 2, there exists an index v (k) for which for all v > Vo(k) such 
that |5j and ^rrjr are integers, there exists a Sterner equiangular tight frame of v{l + |5y) 

1 v(v— 1) 

fc(fc-l) ■ 
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vectors for a space of dimension — — - 



M 


A T 


7 

k 


V 


r 


R/C Construction of the Steiner system 




6 


16 


2 


4 


3 


TTT\ (~\ 111 f* 

R 2-blocks of v — 


4; Affine with g = 2, n = 2 




7 


28 


3 


7 


3 


TTT> 1~\ 1 1 1 f 

R 3-blocks of v = 


7; Projective with q = 2, n = 


= 2 


28 


64 


2 


8 


7 


R 2-blocks of v = 


8; Affine with q = 2, n = 3 




•35 


120 


3 


15 


7 


R 3-blocks ofv = 


15; Projective with q = 2, n 


= 3 


66 


144 


2 


12 


11 


R 2-blocks ofv = 


12 




99 


540 


5 


45 


11 


R 5-blocks of v — 


45 




3 


9 


2 


3 


2 


C 2-blocks of v — 


3 




10 


25 


2 


5 


4 


C 2-blocks of v = 


5 




12 


45 


3 


9 


4 


C 3-blocks of v = 


9; Affine with q = 3, n = 2 




13 


65 


4 


13 


4 


C 4-blocks of v — 


13; Projective with q = 3, n 


= 2 


15 


36 


2 


6 


5 


C 2-blocks of v — 


6 




20 


96 


4 


16 


5 


C 4-blocks of v = 


16; Affine with q = 4, n = 2 




21 


49 


2 


7 


6 


C 2-blocks of v = 


7 




21 


126 


5 


21 


5 


C 5-blocks of v — 


21; Projective with g = 4, n 


= 2 


26 


91 


3 


13 


6 


C 3-blocks of v = 


13 




30 


175 


5 


25 


6 


C 5-blocks of v = 


25; Affine with q = 5, n = 2 




31 


217 


6 


31 


6 


C Projective with q = 5, n = 2 




36 


81 


2 


9 


8 


C 2-blocks of v = 


9 




45 


100 


2 


10 


9 


C 2-blocks of v = 


10 




50 


225 


4 


25 


8 


C 4-blocks of v — 


25 




55 


121 


2 


11 


10 


C 2-blocks of v — 


11 




56 


441 


7 


49 


8 


C Affine with q = 


7, n = 2 




57 


190 


3 


19 


9 


C 3-blocks of v — 


19 




57 


513 


8 


57 


8 


C Projective with q = 7, n = 2 




63 


280 


4 


28 


9 


C Unital with q = 


= 3; Denniston with r = 2, s - 


= 3 


70 


231 


3 


21 


10 


C 3-blocks of v = 


21 




72 


640 


8 


64 


9 


C Affine with q = 


8, n = 2 




73 


730 


9 


73 


9 


C Projective with q = 8, n = 2 




78 


169 


2 


13 


12 


C 2-blocks of v — 


13 




82 


451 


5 


41 


19 


C 5-blocks of v — 


41 




90 


891 


9 


81 


10 


C Affine with q = 


9, ra = 2 




91 


196 


2 


14 


13 


C 2-blocks of v — 


14 




91 


1001 


10 


91 


10 


C Projective with q = 9, n = 2 




100 


325 


3 


25 


12 


C 3-blocks of v = 


25 





Table 2: The ETFs of dimension 100 or less that can be constructed by applying Theorem Q] to the eight 
infinite families of Steiner systems detailed in Section[3] That is, these ETFs represent the first few examples 
of the general constructions summarized in Table [T] For each ETF, we give the dimension M of the 
underlying space, the number of frame vectors N, as well as the number k of elements that lie in any block 
of a u-element set in the corresponding (2, k, u)-Steiner system. We further give the value r of the number of 
blocks that contain a given point; by Theorem[2j |(/ n , f n ')\ = h measures the angle between any two frame 
elements. We also indicate whether the given frame is real or complex, and the method(s) of constructing 
the corresponding Steiner system. 
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In particular, for any fixed k > 2, letting v be either jkik — 1) + 1 or jk(k — 1) + k for 
increasingly large values of j results in a sequence of Steiner equiangular tight frames whose 
redundancy is asymptotically k; these frames can be real if there exist real Hadamard matrices 
of sizes jk + 1 or jk + 2, respectively. 

Proof. To prove the necessary conditions on M and N, recall that Steiner ETFs, namely 
those ETFs produced by Theorem [T], have N = t>(l + an d M = |fez^ ■ Together, these 
two equations imply N = v + kM. Solving for k, and substituting the resulting expression 
into N = v(l + §5±) yields the quadratic equation = (M -l)v 2 + 2(N - M)v- N(N - M). 
With some algebra, the only positive root of this equation can be found to be v = as 
claimed. Substituting this expression for v into N = v + kM yields k = M( ^ +a ^ ■ Having v 
and k, the previously-mentioned relations bk = vr and v — I — r(k — 1) imply r = |5j = 
and b = |r = M, as claimed. 

The second set of conclusions is the result of applying Theorem [T] to the aforementioned 
(2, k, v )-Steiner ETFs that are guaranteed to exist for all sufficiently large v, provided rEr 
and jjfezi) are integers. The final set of conclusions are then obtained by applying this 
fact in the special cases where v is either jk(k — 1) + 1 or jk(k — 1) + k. In particular, if 
v = jk(k — 1) + 1 then |5y = jk and M = jjjfezH = j[jk(k — 1) + 1] are integers, and the 
resulting ETF of (jk + l)[jk{k — 1) + 1] vectors has a redundancy of k + \ that tends to k 
for large j; such an ETF can be real if there exists a real Hadamard matrix of size jk+ 1. 
Meanwhile, if v = jk(k - 1) + k then fry = jk + 1 and M = = (j^ + l)b'( fc -!) + !] 

are integers, and the resulting ETF of fc(j'A; + 2)[j(k — 1) + 1] vectors has a redundancy of 
^jf+f ^ na ^ tends to k for large j; such an ETF can be real if there exists a real Hadamard 
matrix of size jk + 2. □ 

We conclude this section with a few thoughts on Theorems [T] and [2J First, we empha- 
size that the method of Theorem [T] is a method for constructing some ETFs, and by no 
means constructs them all. Indeed, as noted above, the redundancy of Steiner ETFs is 
always strictly greater than 2; while some of those ETFs with |j < 2 will be the Naimark 
complements of Steiner ETFs, one must admit that the Steiner method contributes little 
towards the understanding of those ETFs with = 2, such as those arising from Paley 



graphs |30[. Moreover, Theorem |5] implies that not even every ETF with ^ > 2 arises from 



a Steiner system: though there exists an ETF of 76-elements in M 19 30] , the corresponding 



parameters of the design would be v = y, r = 5 and k = y, not all of which are integers. 

That said, the method of Theorem [1] is truly significant: comparing Table [2] with a 
comprehensive list of all real ETFs of dimension 50 or less [30J, we see the Steiner method 
produces 4 of the 17 ETFs that have redundancy greater than 2, namely 6 x 16, 7 x 28, 
28 x 64 and 35 x 120 ETFs. Interestingly, an additional 4 of these 17 ETFs can also be 
produced by the Steiner method, but only in complex form, namely those of 15 x 36, 20 x 96, 
21 x 126 and 45 x 100 dimensions; it is unknown whether this is the result of a deficit in 
our analysis or the true non-existence of real-valued Steiner-based constructions of these 
sizes. The plot further thickens when one realizes that an additional 2 of these 17 real ETFs 
satisfy the necessary conditions of Theorem [2J but that the corresponding (2, k, v )-Steiner 
systems are known to not exist: if a 28 x 288 ETF was to result of Theorem [lj the 
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corresponding Steiner system would have k = 6 and v = 36, while the 43 x 344 ETF would 
have k = 7 and v = 43; in fact, (2, 6, 36)- and (2, 7, 43)-Steiner systems cannot exist With 
our limited knowledge of the rich literature on Steiner systems, we were unable to resolve 
the existence of two remaining candidates: 23 x 276 and 46 x 736 ETFs could potentially 
arise from (2, 10,46)- and (2, 14, 92)-Steiner systems, respectively, provided they exist. 

4. Equiangular tight frames and the restricted isometry property 

In the previous section, we used Theorem [1] to construct many examples of Steiner ETFs. 
In this section, we investigate the feasibility of using such frames for compressed sensing 
applications. Here, we identify a frame F = {f n }n=i i n with its synthesis operator 
F : C N — > Hjv/, Fg := Yln=i 9( n )fn- That is, F is an M x N matrix whose columns are the 
/ n 's. For a given 5 and if, such an operator F is said to have the (if, 5) -restricted isometry 
property ((if , <5)-RIP) if: 

(l-«< \\Fg\\l < (1 + 8)\\g\\l (1) 

for all g G C N that are K -sparse, that is, for which g(n) ^ for at most if values of n. 
The central problem of compressed sensing is to efficiently solve the underdetermined linear 
system Fg = f for g, given that / itself arises as / = Fg where g is if-sparse. Here, 
the true challenge is that, despite the fact that / is a linear combination of at most if of 
the / n 's, one does not know a priori which particular if vectors were employed. Moreover, 
since the values for N and if encountered in applications are typically very large, it is not 
computationally feasible to check every if-subset of {f n }n=i- ^ * s therefore a remarkable 
fact 0] that g = go can indeed be efficiently recovered from the system Fg = f using linear 
programming, provided the operator F is (2K, 5) -KIP for some 

Since RIP is such an exceedingly nice property, it is natural to ask whether such matrices 
even exist — they do. In fact, [8( used concentration-of-measure arguments to show that 
for every 5 > 0, there exists a constant C such that M x N matrices of Gaussian or 
Bernoulli (±1) entries are (K,5)-K1P with high probability provided M > CK\og(N/M). 
Similarly, matrices formed by taking random rows of a Fourier matrix satisfy RIP with high 
probability when M > CK log 4 (iV) jiij. These existence results have spurred a great deal 
of interest in deterministic RIP matrix constructions that have M = O (if log' 3 N) for some 
> 1, but no such constructions are known to date. Instead, the best known deterministic 
constructions, such as the one given in 0, have M = 0(K 2 ). Despite the fact that all 
ETFs do indeed match this state-of-the-art level of performance, sadly some ETFs — Steiner 
ETFs in particular — fail to do any better. 

To clarify, let K. be any if -element subsequence of {1, . . . , N}, and let Fjc := {f n }neic 
be the corresponding M x K submatrix of F. Using a standard argument, one may show 
that the (if, <5)-RIP condition is equivalent to having the spectrum of each sub-Gramian 
FfcF K lie in the interval [1 — 5, 1 + 5}; in frame parlance, this implies that each Fjc is a 
good Riesz basis. Letting p(A) denote the spectral radius of a given if x if matrix A, the 
(if, 5)-RIP condition is equivalent to having p(F^F K — I) < 5 for all JC. The problem of 
constructing RIP matrices thus reduces to one of spectral estimation. At this point, most 
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research on constructing RIP matrices falls back on a simple, but effective tool: Gershgorin 
circles, namely the fact that every eigenvalue of A lies in one of the disks in the complex 
plane centered at a,k,k and having radius Ylk'^k \ a k,k'\- When the / n 's have unit norm, as 
in the case of ETFs, the diagonal entries of the sub-Gramian F^F^ are all one, and so the 
application of Gershgorin's estimate to A = F^F^ — I reduces to the so-called worst-case- 
coherence bound: 

max P( F K F K ~ l ) = max max V | (f n , f n >)\ < (K — 1) max | (f n , f n> ) \ . 

n'GK 

In particular, in order for F = {f n }n=i to be (K, <5)-RIP, it suffices to have: 

(K-l)max\(f n ,f n ,)\<6. (2) 

Further note that when using to demonstrate RIP for any fixed 5 < 1, the largest possible 
values of K occur when \(f n ,fri)\ achieves its lower Welch bound ( m^v^i) ) 5 ' which occurs 
precisely when F is an ETF |27|. We summarize the preceding discussion as follows: 

Theorem 3. For any fixed 5 < 1, an equiangular tight frame F = {f n }n=i ^ n % has the 
(K, 5) -restricted isometry property ([I]) for all K < 1 + 5( M ^l^ )^ ■ 

Moreover, for any unit norm f n 's, no argument that relies on the worst-case-coherence-based 
bound ([2]) can provide a better range for such K . 

Note that when iV > 2M, we have 1 < — 2 and therefore the maximum permissible 
value of K is on the order of Ms, which is consistent with other known deterministic 
constructions of RIP matrices. This is not to say that ETFs, in general, cannot be RIP with 
M = 0(K\ogP N) for some /3 > 1, but rather, that such a fact cannot be obtained using 
the worst-case-coherence-based bound (|2J). 

This hope notwithstanding, one of the sad consequences of the Steiner construction 
method of Theorem [1] is that we, for the first time, know that there is a large class of ETFs 
for which the seemingly coarse estimate 02]) is, in fact, accurate. In particular, recall from 
Theorem [TJ that every Steiner ETF is built from carefully overlapping v regular simplices, 
each consisting of r + 1 vectors in a r-dimensional subspace of ^-dimensional space. In 
particular, letting K — r + 1, the corresponding subcollection of all K vectors that lie in a 
given block are linearly dependent, which, in accordance with (TjQ), forces the corresponding 
5 to be at least 1. Recalling the value of r given in Theorem [2j we see that Steiner ETFs 
cannot be (K, <5)-RIP for any 5 < 1 so long as K is at least 1 + ( n^~m ) 5 • That is, for 
Steiner ETFs, the best one can truly do is, in fact, given by Theorem |3J This begs the open 
question: are there any ETFs which are RIP with M = 0(K log' 3 N), or does optimizing 
a coarse bound — worst-case-coherence — always come at the cost of being able to realize a 
truly small spectral radius? 
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